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(57) ABSTRACT 

Mobile system and method for monitoring environmental 
parameters involved in growth ormetabolic transformation of 
algae in a liquid. Each of one or more mobile apparati, sus- 
pended or partly or wholly submerged in the liquid, includes 
at least first and second environmental sensors that sense and 
transmit distinct first and second environmental, growth or 
transformation parameter values, such as liquid temperature, 
temperature of gas adjacent to and above the exposed surface, 
liquid pH, liquid salinity, liquid turbidity, 0 2 dissolved in the 
liquid. CO, contained in the liquid, oxidization and reduction 
potential of the liquid, nutrient concentrations in the liquid, 
nitrate concentration in the liquid, ammonium concentration 
in the liquid, bicarbonate concentration in the liquid, phos- 
phate concentration in the liquid, light intensity at the liquid 
surface, electrical conductivity of the liquid, and a parameter 
a(alga) associated with growth stage of the alga, using PAM 
fluorometry or other suitable parameter measurements. 

31 Claims, 5 Drawing Sheets 
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ENVIRONMENTAL MONITORING OF 
MICROBE METABOLIC TRANSFORMATION 

ORIGIN OF THE INVENTION 

5 

The invention described herein was made in the perfor- 
mance of work under a NASA contract and by an employee of 
the United States Government and is subject to the provisions 
of Section 305 of the National Aeronautics and Space Act of 
1958, as amended, Public Law 85-568 (72 Stat. 435, 42 10 
U.S.C. §2457), and may be manufactured and used by or for 
the Government for governmental purposes without the pay- 
ment of any royalties thereon or therefore. 

FIELD OF THE INVENTION 15 

This invention relates to monitoring of a plurality of envi- 
ronmental and growth parameters associated with growth of 
algae and other members of the archaea, bacteria and/or 
eukarya domains (collectively referred to herein as “alga” or 20 
“algae” for convenient reference). 

BACKGROUND OF THE INVENTION 

Growth or enviromnental control or other controlled meta- 25 
bolic transformation of archaea, bacteria and/or eukarya 
domains in a liquid environment is challenging, in part 
because the ambient liquid allows transport of nutrients and 
other substances between the microbes, whether of the same 
species or of different species. Monitoring of the individual 30 
microbes for one, two or several changes in environmental, 
growth and/or metabolic transformation parameters is often 
required, if optimal control is desired. 

Monitoring and control of a large number of environmental 
and growth parameters is often required in order to respond to 35 
predictable and unpredictable events that occur during 
growth, transformation and harvesting of algae. These activi- 
ties should permit distinction between parameter values at 
different locations and different times in the growth and meta- 
bolic transformation processes, should permit optional pre- 40 
scription of vehicle trajectories, should provide for collection 
and prompt analysis of the measurements, and should pemiit 
use of multiple vehicles. 

SUMMARY OF THE INVENTION 45 

These needs are met by the invention, which provides a 
method of monitoring selected environmental and/or growth 
parameters for growth of algae and other archaea, bacteria 
and/or eukarya immersed in a liquid environment. The 50 
method includes the following processes: (1) one or more 
distinct organisms is provided in a liquid; (2) a mobile system 
that floats on and/or is submergible in the liquid is provided, 
the mobile mechanism comprising at least first and second 
environmental sensors that sense distinct first and second 55 
environmental parameter values, pv(t;l) and pv(t;2), respec- 
tively, that are drawn from liquid temperature, temperature of 
gas adjacent to and above the exposed liquid surface, liquid 
pH, liquid salinity, liquid turbidity, 0 2 dissolved in the liquid, 
C0 2 contained in the liquid, oxidation/reduction potential of 60 
the liquid, nitrate concentration in the liquid, ammonium 
concentration in the liquid, bicarbonate concentration in the 
liquid, phosphate concentration in the liquid, other nutrient 
concentrations in the liquid, light intensity at the liquid sur- 
face and below the surface, electrical conductivity of the 65 
liquid, and a parameter r|(alga) associated with growth stage 
of the alga; (3) the first value pv(t;l) is sensed at spaced apart 
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first and third sensing times, t=tl and t=t3; (4) the second 
value pv(t;2) is sensed at a second sensing time, t=t2, that is 
intermediate between the first and third sensing times ; and (5 ) 
the sensed values, pv(tl;l) and pv(t3 ;1), sensed at the first and 
third sensing times, and the value pv (t2;2), sensed at the 
second sensing time, are transmitted to a receiver that is 
spaced apart from the mobile system. Pulse amplitude modu- 
lation (PAM) fluorometry or another suitable mechanism is 
applied for some of these measurements. 

Where a plurality of mobile systems are provided, the 
sensed values of a selected parameter at different locations 
are optionally combined with each other and a function, 
dependent upon location and time and representing a field of 
values of this parameter within or adjacent to the growth 
medium, is constructed for one or more sampling times and 
optionally displayed for a user. Optionally, this field of values 
is analyzed to estimate relevant enviromnental and/or growth 
parameter values and trends, to compare at least one such 
value with a range of acceptable values for such parameter, 
and to determine whether any change(s) should be introduced 
for such parameter. 

Optionally, the mobile system is allowed to move with 
prevailing liquid currents in the liquid, and measurements of 
one or more parameter values are sampled and transmitted at 
each of a sequence of spaced apart times and locations, 
together with an estimate of the sampling location coordi- 
nates, angular orientations and corresponding sampling 
times. Alternatively, the mobile system is provided with a 
location sensing and/or angular orientation sensing mecha- 
nism and with a locomotion mechanism that is configured (i) 
to move the mobile system along a selected trajectory in the 
liquid and/or (ii) to alter or correct the mobile system present 
location to a selected corrected location and/or (iii) to alter or 
correct the mobile system present angular orientation to a 
selected corrected orientation. 

Optionally, a time rate of change of at least one parameter 
value pv(t;l) is estimated, based on an elapsed time differ- 
ence, At=t2-tl, and this elapsed time difference is adjusted 
subsequently according to how quickly or slowly the param- 
eter value pv(t;l) is changing with time. 

Optionally, the mobile system is provided with a submer- 
sion mechanism that permits the mobile system, or a portion 
of the system containing at least one of the sensors, to be 
submerged by a selected, continuously or discretely variable 
depth within the liquid. 

Optionally, the mobile system has a structure with at least 
one aperture, through which at least one sensor performs 
measurements of an associated enviromnental parameter. 
Optionally, measurements of two or more parameters are 
performed through this aperture. Optionally, at least one 
sample of liquid, algae or other relevant substance can be 
captured through this or another aperture and preserved for 
subsequent analysis. 

On-board, or subsequent, analyses can be implemented 
using procedures such as PAM fluorometry of the samples to 
determine photosynthetic efficiency or other parameters. 

Optionally, where N mobile systems are provided (N>2), a 
mutual repulsion system is provided for each mobile system 
so that any two such mobile systems that come within a 
selected distance of each other will repel each other, in order 
to prevent “bunching” or congregation of two or more such 
mobile systems within a small region. 

The invention described in the preceding can be applied to 
environmental, growth and/or metabolic transformation 
monitoring of algae, in a liquid environment. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 illustrates use of the invention, with one mobile 
system, in a liquid environment. 
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FIGS. 2 and 3 illustrate an embodiment of a submersion 
and flotation mechanism used in the invention. 

FIG. 4 illustrates a mutual repulsion system for two or 
more mobile systems. 

FIG. 5 illustrates construction of a field of values deter- 
mined according to the invention. 

FIG. 6 illustrates estimation of location/angular orientation 
coordinates. 

DESCRIPTION OF BEST MODES OF THE 
INVENTION 

FIG. 1 illustrates an environment where the invention can 
be used. A mobile system 11 floats or is submerged in a liquid 
12 having a liquid surface 12LS (exposed or defined by a 
container surface, such as a pipe) and is intended to monitor 
two or more environmental, growth and/or metabolic trans- 
formation parameters (e.g., population size, population com- 
position or metabolic activity) of the liquid or of a substance 
S suspended or otherwise contained in the liquid. The sub- 
stance S may be one or more alga or another member of the 
archaea, bacteria or eukarya domains that is to be grown or 
otherwise transformed over a time interval. Two or more 
environmental and/or growth and/or metabolic transforma- 
tion parameters of the substance S are to be monitored by the 
mobile system, as a function of time and/or as a function of 
location in or adjacent to the liquid 12. These monitored 
parameters may include, but are not limited to, liquid tem- 
perature, temperature of gas adjacent to and above the liquid 
surface, liquid temperature at a selected depth below the 
surface, liquid pH, liquid salinity, liquid turbidity, O, dis- 
solved in the liquid, C0 2 contained in the liquid, oxidization/ 
reduction potential of the liquid, concentration of ammonia, 
nitrate, phosphate, carbonate or other nutrient(s) or inorganic 
substances in the liquid, light intensity where light enters the 
liquid, light intensity at one or more selected depths within 
the liquid, electrical conductivity of the liquid, and one or 
more parameters a(alga) associated with growth stage of the 
alga. The mobile system 11 may float on or be partially or 
fully submerged in the liquid or may be located within and/or 
move within a pipe or other container for the liquid, or the 
mobile system may float on or be submerged relative to an 
exposed liquid surface. 

A first parameter value pv(t;l) is measured at two or more 
spaced apart times, t=tl and t=t2 (>tl), and a second param- 
eter value pv(t;2) is measured at one or more third times, t=t3, 
with tlSt3<t2, and the measured values are transmitted to a 
measured signal receiver 13 that is spaced apart from the 
mobile system 11, using a measured signal transmitter 14 
located on the mobile system. Where M (=2) environmental 
parameters are monitored, times of transmission of the M 
parameters may be interspersed with each other, optionally in 
a round robin format, so that a measured signal representing 
each parameter value is received at one or more spaced apart 
times. Optionally, one or more of the measurements may use 
fluorometry or a similar mechanism. 

Each mobile system 11 is floatable on and/or submergible 
in the liquid 12 or may be located within and/or move within 
a pipe or other closed environment containing the liquid. The 
mobile system 11 optionally includes a submersion mecha- 
nism 31, indicated in more detail in FIG. 3, that permits a 
selected portion of (or all of) the mobile system to be sub- 
merged relative to a liquid surface 12LS to a selected, con- 
tinuously variable or discretely variable depth within the liq- 
uid 12. A submergible portion of the mobile system 11 
optionally includes one or more parameter value sensors 
23 -m (FIG. 1; m=l, 2, . . . , M; Mg2) that each sense one or 
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more relevant environmental, growth and/or metabolic trans- 
formation parameters of the liquid 12 and/or of the substance 
S in the liquid. The measured sensor values pv(t;m) (m=l, 
2, , M) are transmitted by the transmitter 14 to the receiver 
5 13 for analysis and/or storage. Optionally, each sensor 23-m is 
located in or adjacent to one or more apertures 24 -m in a 
containment structure 25 of the mobile system 11, through 
which the sensor measurement can be performed. Optionally, 
two or more sensors are located in or associated with mea- 
111 surement through a common aperture. 

Optionally, each mobile system 11 has a location mecha- 
nism 61 (FIG. 6) that senses the present location coordinates 
(x,y,z), angular orientation (0;4>) (optional), and present time 
15 t of a designated component of the mobile system and pre- 
sents this information to the transmitter 1 4 for transmission to 
the receiver 13 . Optionally, a supplementary location mecha- 
nism is provided to measure or estimate a depth d(m) of the 
sensor 23 -h; below the liquid surface 12LS, for transmission 
20 to the data receiver 13. 

A submersion mechanism 28 -m. which may be but need 
not be the same as the mobile system submersion mechanism 
31, is optionally provided for separate depth control of the 
sensor 23->». The submersion mechanism 28 -m allows the 
25 mobile system 11 to vary the depth of the associated sensor 
23 -in to facilitate sensor readings at arbitrarily selected 
depths d(m) within the liquid 12. 

One implementation of a submersion mechanism 31 pro- 
vides K inflatable bags, 32-A'(k=l,2, . . . ,K;K>2) or other gas 
30 enclosures at spaced apart locations on the mobile system 11, 
or on the sensor 23 -n i , The mass of the mobile system 1 1 , with 
or without one or more selected masses attached, is assumed 
to be chosen so that, when each of the enclosures 32 -k 
includes a selected amount of gas (e.g., air) and has an asso- 
35 ciated baseline volume V0(£), the mobile system and submer- 
sion mechanism 31 floats on the liquid surface 12LS of the 
liquid, in accord with Archimedes principle. When a selected 
amount of gas is released from a one, two or more of the 
enclosures (e.g., 32-1 and 32-2), the enclosure gas volumes 
40 are reduced to V1(1)<V0(1) and V1(2)<V0(2). As a result of 
this volume reduction, the mobile system 11 and/or one or 
more sensors 23 -m are (further) submerged to a controllable 
depth dl below the liquid surface 12LS of the liquid, where 
the depth dl depends upon the volume differences, V0(1)-V1 
45 ( 1 ) and V0(2)-V1(2). 

When the mobile system 11 and/or a sensor 23 -m are to be 
brought closer to the liquid surface 12LS so that the depth d of 
submersion is reduced to d2<dl, gas from a re-pressurization 
mechanism 33 is (re)admitted into one or more of the enclo- 
50 sures 32 -k (not necessarily the enclosures, 32-1 and 32-2, that 
were originally partly deflated), with corresponding increases 
in the enclosure volumes, and the mobile system 1 1 and/or the 
sensor 23 -in rise to a new selected depth d2 (<dl). Appendix 
A contains an analysis of the dependencies of submersion 
55 depth d on enclosure volume differences. 

Optionally, one or more of the sensors 23 -m is located 
adjacent to, or within, an optically or electromagnetically 
transparent aperture 29 -m, through which the sensor 23 -in can 
measure or otherwise sense the environmental parameter(s) 
60 pv(t;m) associated with this sensor. Each parameter value 
measured or sensed by the sensor 23 -m is presented to the 
transmitter 14, for possible transmission to the data receiver 
13, together with a time t(meas) at which the parameter was 
measured or otherwise sensed. Optionally, the measured or 
65 sensed parameter value and associated measurement time 
pv(tmeas;m);t(meas;m)) are transmitted to the data receiver 
13 only if (1) a measurement time difference satisfies 
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?2(meas);m-U (meas;m)^a threshold time difference 

A t(thr;m), (1) 

which may depend upon m and/or (2) a parameter value 
difference satisfies 

\pv{t2;m)-pv{t\;m) I i^a threshold pv difference pv(thr; 

m), ( 2 ) 

which may also depend upon m. 

Preferably, the parameter value and time value signals pv(t- 
meas;m),t(neas;m)) are transmitted using a low power, low 
data rate protocol, such as Zigbee, whose protocols are partly 
specified by I.E.E.E. 802.15.4 standard. Zigbee uses radio 
frequency signals and allows relatively long battery life and 
secure networking. Zigbee operates in the industrial, scien- 
tific and medical (ISM) radio bands, circa 868 MHz in 
Europe, circa 915 MHz intheU.S. and Australia, and circa 2.4 
GHz in most jurisdictions throughout the world. Zigbee can 
be activated in 15 msec or less so that an transmission of a 
isolated signal, such as pv(t(meas);m), can be preceded by 
transmission activation and followed by transmission de-ac- 
tivation, if desired, with resulting low latency and low average 
power consumption. Typical Zigbee applications include 
commercial building monitoring and control, home aware- 
ness and security, home entertainment and control, industrial 
plant management and control, and mobile services, such as 
health care/monitoring and tele-assistance. 

A Zigbee communication system requires a Zigbee coor- 
dinator, one or more Zigbee routers (optional here) for for- 
warding data to other devices, and a Zigbee (non-forwarding) 
end device. The network can form a single cluster for a single 
mobile system 11, or a mesh of clusters (useful for a plurality 
of mobile systems 11). Beacon-enabled networks and non- 
beacon-enabled networks (number of beacons S15) are sup- 
ported. In a non-beacon-enabled network, a non-slotted code 
division multiple access/collision avoidance (CDMA/CA) 
channel access mechanism is used. A beacon-enabled net- 
work provides wakeup/transmit/sleep time intervals that are, 
in one version, multiples of 15.36 msec for transmission rates 
of 250 kbits/sec, to multiples of 24 msec for 40 kbits/sec, to 
multiples of 48 msec for 20 kbits/sec. The I.E.E.E. 802.15.4 
standard specifies the physical layer and the media access 
control layer. Message acknowledgments, when used, do not 
rely upon CDMA/CA. 

Zigbee routers, if any are used, are typically continuously 
active so that a robust power supply is required. The receiver 
(typically located on-shore and non-portable) can be config- 
ured to receive continuously, with the aid of a large(r) battery 
or power supply. The individual mobile systems can transmit 
in a slotted fonnat, with each of M mobile systems typically 
being allotted a slot of length At(mob)«200 msec in a round 
robin cycle of length several sec. Not all environmental 
parameter values pv(t';m) need be transmitted within each 
time slot. Optionally, only the parameter values pv(t;m) that 
have changed significantly since the preceding slot for the 
given mobile system will be transmitted within the present 
time slot, as indicated in Eq. (2). 

Optionally, the housing for the mobile system also includes 
one or more sampling apertures 30 -q (q=l, . . . , Q; Q= 1 ) that 
facilitate physical capture of one or more samples of the 
underwater environment (e.g., a sample of the liquid, of an 
alga, or of surface material in which the alga is grown). The 
sample, after capture, is deposited in and stored in a sample 
repository within the mobile system, for subsequent analysis. 

Where two or more mobile systems, denoted 11-1 and 11-2 
here (FIG. 4) are present, one concern is keeping the mobile 
systems spaced apart from each other by a distance D that is 
at least a minimu m distance D(min) For small but positive 
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values of D(min), for example, D(min)=l-5 cm, each mobile 
system is provided with at least three, preferably at least five, 
mutual repulsion mechanisms, 17 -j-k (j=l, 2; k=l, 2, 3) 
mounted on each such system that provide mutual repulsion, 
5 mobile system to mobile system. One implementation of the 
mutual repulsion mechanisms 17-j-k provides three or more 
magnetic poles, all south poles or all north poles so that any 
two poles, one on mobile system 11-1 and one on mobile 
system 11-2, will repel each other. The three south poles 
to 17-l-A:on the mobile system 11-1 (and, similarly, the three 
south poles 17-2-k on the mobile system 11-2) are preferably 
connected to a single north pole 18^' of an opposite magnetic 
sign, which is optionally located in a region elsewhere within 
the mobile system, as illustrated in FIG. 4. Choice of the 
15 magnetic strength of the pole 1 8-/ will strongly influence s 
distance of closest approach of two mutually repulsive mag- 
nets (both south poles or both north poles), one located on 
each of the two mobile systems. By spacing P south poles on 
a perimeter at approximately uniform angle differences (e.g., 
20 at Aij)=<360 o /P for P such poles), the mutual repulsion of like 
poles will be approximately equally effective, for approach of 
a first mobile system to a second mobile system from any 
direction. 

Estimation of location coordinates, angular orientation 
25 coordinates and corresponding sampling times (x,y,z;0;(|);t) 
can be implemented using GPS signals, for a relatively large 
body LB of liquid, or may be implemented using three or 
more pseudolite signal generators, located at spaced apart 
positions adjacent to the boundary LB, as discussed in detail 
30 in Appendix B. 

Appendix A. Mechanism For Control of Submersion Depth. 

Consider first and second bags or other enclosures 32-1 and 
32-2 (FIG. 3) that is inflatable with a selected gas, preferably 
having a very low density relative to the density of the liquid 
35 12. The enclosures, 32-1 and 32-2, are attached to a mobile 
system 11 and are attached to a controllable, pressurized 
source 33 of the selected gas, which can be introduced into 
one or both of enclosures or can be withdrawn from one or 
both of enclosures. 

40 FIG. 2 illustrates a cross section of a bag or other enclosure 
E, which is assumed to be symmetric about a horizontal plane 
n (referenced as a coordinate value, z=-z p ) and is assumed to 
be (approximately) rotationally symmetric about a z-axis 
(x=y=0). A cross section of the lower and upper sections of 
45 the enclosure shown in FIG. 2 are assumed to be defined as 


z=-zp±flr). 

(A- 1) 


(A- 2) 


50 where r is a radius measured in the horizontal plane n. The 
function f(r) is assumed to be continuously differentiable so 
that the derivative or slope df/dr is defined and continuous for 
OgrSr(max) (edge of the enclosure in the plane n). Note that 
the corresponding slopes {df/dr} (upper) and {df/dr} (lower) 
5 5 on the upper section ( z>-z p ) and on the lower section (z<-z p ) 
of the enclosure are related by 

{4/2rfr}(upper)=-{(#?tfr}(lower). (A-3) 

The total liquid force on the enclosure, assumed to be 
60 directed normal to the enclosure surface at each point, is given 
by 

Gzp)=pOgf 0 K '”‘“ , {{/z p -/r)](^;-)(upper)+{[z p +/fr)] 
{df/dr){\osNQr)}2nr dr/{ 1 +(df/dr ) 2 } 1/2 ,= 

65 pOgf 0 r(max ^ 2z P j {df/dr) (upper) 2jtr dr, (A-4) 

sin ty=(df/dr)/{ \+{df/dr ) 2 } 1/2 . (A-5) 
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where it is assumed that the horizontal plane II is located at a 
depth Az=z p below the surface of the liquid 12, and that the 
mass or pressure equivalent of this liquid (per unit depth 
below the liquid surface) is pO (e.g., 62.4 lbs. per cubic foot 
for water as the liquid 12), A corresponding volume of the 5 
enclosure interior is 

V=j 0 r(max> 2J(r)2jir dr. (A-6) 

The total liquid force F is directed upward, toward negative 10 
values of the coordinate z. Note that, for any smooth shape f(r) 
of the enclosure cross section, the total liquid force F on the 
enclosure is proportional to the depth parameter z p and is 
usually directed upward. Equation of the force F and the 
mobile system weight w(mobile) will provide an equilibrium 15 
depth z P at which the (submerged) enclosure plus weight will 
settle. As the amount of the selected gas introduced into the 
enclosure E increases, the shape function f(r) changes, the 
volume within the enclosure increases, and the equilibrium 
value of the depth parameter z P decreases so that the mobile 20 
system 11 rises in the liquid; and conversely. The net energy 
required to cause the depth parameter z p to decrease by a 
controllable amount so that the mobile system 11 rises in the 
liquid is substantially equal to the energy required to intro- 
duce the additional gas into the enclosure interior, which is 
relatively small. Optionally, when the depth parameter z P is to 
be increased so that the mobile system sinks by a controllable 
amount, a portion of the selected gas in the enclosure interior 
can be bled off and allowed to pass into the ambient medium, 30 
with little or no additional energy expenditure. 

Assuming that the mass of an attached mobile system, 
including ballast masses, if desired, is m(total), the 
enclosure(s) E will rise in the supporting liquid to an enclo- 
sure depth Az=z p for which 

F(z p )=m(total)*g. (A-7) 

This relation may be corrected, if desired, for the supplemen- 
tal buoyancy for the mobile system 1 1 that is provided by the 
supporting liquid. 

Appendix B. Estimation of Location Coordinates and Angu- 
lar Orientation 

Coordinates and Sampling Times. 

Consider a mobile system 61 -m (m=l, .... M;Mgl) that 45 
is floating on or submerged in a liquid 62, having a boundary 
3B, as illustrated in FIG. 6. Three or more spaced apart 
pseudolites, 63-A(k=l, . . . , K; KS5) are located at selected 
distances from the boundary 3B, or optionally within the 
liquid, each pseudolite 63-Ahas known and fixed coordinates 50 
(x k ,y k ,z k ), and any three of the pseudolites are non-collinear. 

A signal receiver 62 is located at a selected position, relative 
to or on the mobile system 61-??/. The mobile system 61-?// has 
a timing mechanism that is coordinated with a timing mecha- 
nism on each of the pseudolites 63-A. At each of a specified 55 
sequence of spaced apart times, {t=t im )} m , the pseudolite 
63 -A transmits a short location/timing signal S L/7 (t km ), pref- 
erably including a coded portion that identifies the pseudolite 
source and specifying the time this signal is transmitted. The 
signal S^/^t^,) is received by a receiver antenna on the 
mobile system 61-?//, which computes a time difference 

A A- (B-l) 

between a time, t=t Mji , this signal is received by the receiver 
antenna (m) and the signal transmission time for the pseudo- 
lite (k). A separation distance, 


d t „{t)=cAt t . m (B-2) 

= { Mr))-x k (p)f + (>„(?)) - ydp)? + 
te„ W)-/*(P)) 2 1 1/2 , 

between the receiver antenna, having location coordinates 
(x m (r),y m (r)),z m (r)), and the pseudolite, having location coor- 
dinates (x t (p),y t (p)),z t (p)), is then estimated, where c is the 
propagation velocity of an electromagnetic signal in the 
ambient environment (the “speed of light”). 

Each pseudolite will have fixed location coordinates that 
do not vary with time, and preferably any temporal variation 
between the timers foreachpairofpseudolites is immediately 
corrected, through p-to-p communication. Flowever, an indi- 
vidual receiver associated with a mobile system (m) will have 
varying location coordinates, and the on-board timer for this 
mobile system may have a time shift, At=x(t;m) that varies 
with time and varies with the particular mobile system (m) 
being examined. This time shift (positive, zero or negative; 
unknown initially) is explicitly incorporated in the equations 
for distance as follows: 

{(xjr))-x k ip)) 2 + iyjr))-y t ip)) 2 +(z m (r))- 

z k (p)) 2 } ll2 = c {At(k. m )+r(t;m)}(k=l , .... K;m 

fixed), (B-3(k,m)) 

where At(k,m) is a time difference measured at the pseudolite 
number k or measured at the mobile system number m. Equa- 
tions (B-3(k,m)) fork=2, 3, 4, etc. (m fixed) are squared and 
subtracted from Eq. (3-B(k=l pm) as follows. 


1 (Xm(r))-X k (p)f + (y„(r)) - y k (p)) 2 + (z m (r)) -ft (p)f}~ 

{ Mr)) - Xi (p )) 2 + (y„(r)) - yi(p )) 2 + (z m (r)) - zi {p)f} = 
liXmWx^p) - x k (p)) + y„(r)(yi(p) - y k (p)) + z m (r)(Zi(p) - Zk(p))) + 
x k (pf -*i(p) 2 +y k (p ) 2 -y k (p ) 2 +Zk(p ) 2 -Zi(p ) 2 += 
c 2 {At\k.m) + r(r; m)} 2 -c 2 {Ar(l.m) + r(r; m)} 2 += 
c 2 {2r(t; m)(At(k , m) - A/f 1 , mj) + Alik., m) 2 - A/(l, - 4(k. m)) 


The unknowns here are the location coordinates, x m (r), y m (r) 
and z m (r), and the time shift t(t;m) for the mobile system 
61 -m; all other quantities are known. Equations (B-4(k,m)) 
can be re-expressed in matrix form as 


MX = H. 




(B-5) 

xi (p)-x k (p) 

(yi(p)-y k (p) 

zi {p)-z k (p) 

A t(k, m) - " 
Ar(l, m 

(B-6) 

xi(p)-x k (p) 

(yi(p)-y k (p) 

Zl(p)~Z k (p) 

A t(k, m) - 
Af(l, m) 


xi(p)-x k (p) 

(yi(p)-yt(p) 

Zl (P)-ZkiP) 

A t(k, m) - 
Af(l, m) 


Xl(p)~X k (p) 

(yi(p)-y k (p) 

Zl (p)-z k (p) 

A t(k, m) - 
Af(l,m) 
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-continued 

X ,r = [x m (r)y m (r) z m /r) r(r,m)], (B-7) 

2 H" = [ x ! ( p ) 2 - x k (p) 2 + y, (p) 2 - 

»s(+Zi(p) 2 -Zt(p) 2 + c 2 (At(k, m) 2 -A /(/, m) 2 ), x k (p) 2 - 

Jtt(p ) 2 +yi(p ) 2 -yd+zdp ) 2 -zi(p ) 2 + 
c 2 (At(k, m) 2 - A*/, m) 2 ), * (p) 2 - x t (p) 2 + yi (p) 2 - 
>>j(+Zl(P) 2 -Zk(p) 2 Fc 2 (At(k, m) 2 -AtU, m) 2 ), 
xi (p) 2 -x k (p) 2 +yi(p) 2 - 
yk (+Zi (pi 2 - Zi (pi 2 + c 2 (At(k, m) 2 - At (/, m) 2 )] 


erably, KS20, and more preferably KsilO), a three-dimen- 
sional field or surface of values that precisely matches the 
values pv(x t ,y i ,,z (t ,t t ;ni) for the coordinates (x=x k ,y=y k ,z=z k ) 
is constructed, as illustrated in FIG. 5. This field of values is 
5 defined by 


V* K 

F(x,y,z;t*tO)y pv(x k ,y k ,z k ,t k ,m) Y\ {(*- 


(C-l) 


The matrix M has an inverse M -1 (rank 4 ), in part because no 
three of the location coordinate triples (x^pfy^pfz^p)) (i=l , 

2, 3, 4, 5)forthepseudolites are collinear. The matrix solution 
X of Eq. (B-5) becomes 15 

X=M~ l H. (B-9) 

The time shift x(t;m) for mobile system 61 -m may, but need 
not, vary with time t. 

The location coordinates (x„,(t;r),y m (t;r),z m (t;r)) estimated 20 
for mobile system 61 -m will vary with time t as the mobile 
system moves on or within the liquid 62. Two consecutive 
sampling times, t„ and t„ +1 , for the mobile system 61-»/ may 
be determined by identifying a time difference, t„ +1 -t„, for 
which a separation distance first satisfies 25 

{(a„i.xFa,x)) 2 ^,(w'>-^) 2 + 

( z m(t n+l ;r)-z m (t„;r)) 2 } v2 ^d(thr), (B- 1 0) 

where d(thr) is a selected threshold separation distance, 
which is chosen according to an estimated difference in one or , f) 
more environmental parameters that may occur with the sepa- 
ration distance d(tlir). This sampling frequency may vary 
from one environmental parameter to another. 

Angular orientation coordinates (<|>,0)=(azimuthal angle, 
polar angle) for the mobile system, relative to a horizontal , 
plane may be estimated by obtaining the location coordinates, 

(x m (r; bow) ,y m (r; bow) (r;bow)) and 

(x w (r;stem),> , w ( 2 - ;stem)/ m (r;stern)) 

for spaced apart, first and second locations on the bow and on 40 
the stem, respectively, of the mobile system 61-»;. These 
location coordinates may be estimated, applying the consid- 
erations involved in Eqs. (B-l) through (B-9) to measure- 
ments made at a second data receiver located on the bow and 
at a third data receiver located on the stern of the mobile 45 
system 61-m. The angles 4>,„ and 0 OT for mobile system no. nr 
are then estimated from the relations 


Xj )(y - yjHz - zj)/[x t -xjXyj - y k )(zj - z k )}, 
F(x k ,y k ,z k , ; r » (0) = pv(x k , y k ,Zk. t k \m), (C-2) 


where, for each value k, the finite product II is fomied for all 
ratios (x-x )/(x k -x ) for which j^k. The function F is a poly- 
nomial of degree K in each of the coordinates x, y and z and 
is representable as a surface S of height pv(x A .,y Ar ,z Ar ,t i .;m) 
above a plane II at the location (x k ,y k ,z k ), as illustrated in FIG. 
5 . 

Where two or more measurement locations (x k ,y k ,z k ) have 
at least one of the same coordinate values, x k or y k or z k , a new 
“rotated plane” with coordinates (x',y',z') is introduced, 
defined by 


X 


V 

y 

= /W 2 (j8)-M 3 (J 

/ 

. z . 


t! 


Mi(a) = 


cos a sino’ 0 
— sino cos a 0 


0 0 1 


cos/3 0 sin j3 


Mm = 


0 1 0 
-sin /3 0 cos/? 


(C-3) 


(C-4) 


(C-5) 


The angular variables a and |) are chosen so that an equivalent 
finite product 


* K 

7, (x',y',z/;txtO) \ pv(x' k ,y' k ,z' k ,t k \m) ]”| ((*' - 

4-H j=i,M 


(C-6) 


-y'j)(z. -z))/(4 - x'jKy'j - y' k )Zj - z k ) 1, 


tan(p m = {^ OT (r; bow) — y m (r; stern)} / (/"; bow) — x m (r; stem)}, (g_ j 

tand m = { z m (r ; bow) - z m ( r 'i stem)} /{x m (r, bow) - x m (r; stem)) 2 + (B- 12) 

(y„(r; bow) - y„(r; stem)) 2 ) 1 ' 2 , 

Nonnally, the polar angle 0 m will be close to 0, because the 
mobile system will be approximately horizontal in the liquid 
62. 

Appendix C. Construction of a Surface From a Field of Val- 
ues. 

Consider a field of measured values ^y{x k ,y k ,z k ,X k \m) 
(k=l , . . . , K) of a selected environmental parameter (m fixed), 
where the values have been measured at times (t^—tO, a time 
value independent of k) that are reasonably close to each 
other. This is preferably done using two or more mobile 
systems, which may be individually stationary or may move 
relative to each other. If the integer K is not too large (pref- 


50 

expressed in terms of the equivalent coordinates (x',y',z'), 
does not include any situation in which x'j=x' k or y'j=y\ or 
z'=z' k . This can be achieved in a straightforward maimer, 
because the total number of coordinates (x 7 ,y p zJ) is finite and 
55 equal to 3K. The coordinates (x',y',z') are then expressed in 
terms of the original coordinates (x,y,z) as 


x' 


X 

y 

= M 3 (-ar)M 2 (- 1 S) 

y 

z' _ 


. z . 


which is the inverse of Eq. (C-3). The resulting field of values 
65 function F"(x,y,z;t=t0) will be a polynomial with the same 
maximum degree (K) in each of the coordinates but may not 
be expressible in the elegant form set forth in Eq. (C-l). In 
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order to avoid high polynomial degrees, with the associated 
frequency variation in each of the location coordinates, it is 
preferable to limit the number K of sampled pv values to a 
relatively small integer, such as K=3-6. 

It is important to control certain environmental parameters 
to insure that these parameters stay within a suitable range in 
order to avoid adversely interfering with the growth or meta- 
bolic transformation process for a particular microbe. These 
parameters include: liquid pH (preferred range=3.5SpH 
(liq)S9.5); liquid temperature T (0° C.£=T(liq)S40° C.); 
C0 2 concentration (0-10 mmoles/liter); 0 2 concentration 
(0-1 mmole/liter); nutrient concentration (0-1 mmole/liter); 
light irradiance I adjacent to a liquid surface (0-2500 photons/ 
(meter)-2-sec). 

Appendix D. Application of Fluorometry to Transformation 
Measurements. 

In a fluorescence process, incident light having energy 
E=hv=hc/X is absorbed by an object, a portion of the light 
energy is transformed into one or more other processes, and 
the remainder, if any, is re-emitted at a higher wavelength, 
k’>k or v'<v, representing a lower emission energy. The 
remainder of the incident energy may be absorbed and used 
for photosynthesis, for heat dissipation and/or for quantum 
processes such as (temporary) atomic or molecular excita- 
tion. This perspective may be expressed as an equation 

£=hv+photosynthesil (£ /?s )+heat dissipation(E <Ass )+ 

quantum processes(£‘^ p )+fluorescence ( E FL =hv ') (D-l) 

For simplicity, the presence of quantum processes is ignored 
here. 

When a leaf or other biological object is kept in the dark, 
with no energy supplied, the characteristic internal energy of 
the object will decrease toward a minimum energy value E mjn . 
When this object is illuminated with incident light (E hv), the 
reissued fluorescence signal energy E FL will increase to a 
maximum value E max , with the remainder, E ps +E diss , being 
generally not observable optically. The difference, E,,, = 
E max -E min , is referred to as the variable fluorescence and has 
a range of values, E v «0.8 E max , for healthy, growing plants. 
For many algae, the ratio EJE max has an optimal range, 
depending upon the process of interest., and the system 
approaches this optimal range monotonically as growth or 
other metabolic transformation proceeds. 

When non- saturating light is received by a target (e.g., 
chlorophyll a), the target is elevated to an excited state in 
which an electron is transferred into an electron transport 
chain (ETC), which can be subsequently used for other pur- 
poses, such as production of ATP, reduction of NADP, and 
production of glucose. After electron transfer to and through 
an ETC has occurred, the electron acceptor becomes avail- 
able again to accept another electron. When higher light 
intensity is delivered, the reaction centers become tempo- 
rarily inactivated (closed) until electrons presently in the 
reaction centers are cleared out through transfer. When this 
saturation and temporary closure occurs, excess incident light 
is re-emitted as fluorescence, which can be measured and 
used as a diagnostic parameter to estimate the present status 
of the target. When a group of cells is under stress, saturation 
occurs more easily, and the system processes the incident 
light less efficiently than would occur in a group of cells not 
under such stress. 

A fluorescence target is initially kept in the dark, to allow 
the ETC to become cleared of electrons and to open all accep- 
tor channels. At this point, the system is most efficient in 
accepting small amounts of incident light and in processing 
the electrons produced in the excited states. As the system 
approaches saturation, efficiency decreases. System effi- 
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ciency changes more quickly when system stress is present 
and saturation occurs at reduced (time-integrated) incident 
light levels. 

Photo-excitation of the fluorescence target is preferably 
5 implemented using pulse light, with an excitation interval 
having a temporal length in a range of 0.1 psec-1 msec, in 
order to avoid saturating the target and to estimate the mini- 
mum value Delivery of a single, longer pulse, with an 
associated time interval length in a range 0.1-2 sec can be 
111 used to estimate a fluorescence saturation eneigy level. 

A pulse amplitude modulation (PAM) fluorometer can 
deliver incident light pulses in one or several different wave- 
length ranges (X) and at different light delivery rates r(X), in 
1 . order to estimate the efficiency parameter F y/F M (comparable 
to EJE max ) and an initial value E min . The fluorescence effi- 
ciency will vary with the wavelength range used and with the 
light delivery rate. A PAM fluorometer uses three or four 
wavelength ranges to estimate the presence and concentra- 
20 tions of at least three groups of algae: (i) green algae (sensed 
by presence of chlorophyll b); (ii) cyanobacteria (sensed by 
presence of allophycocyanin), and (iii) diatom/dinoflagellate 
(sensed by presence of chlorophyll c, fucoxanthin and caro- 
tenoids). Other wavelength-related groups can also be used 
25 here. Stress of a group of cells in the target can be imple- 
mented through O, deprivation, CO, deprivation, nitrate 
(salt) excess or deficit, or other nutrient excess or deficit, 
among others. 

PAM fluorometry allows control of at least four parameters 
30 for monitoring a present stage of micro-organism and 
selected environmental parameters: wavelength range, light 
intensity or irradiance, duration of exposure, and duty cycle 
(fraction of time the micro-organism is exposed to light). 

35 What is claimed is: 

1 . A method of monitoring an environment for growth or 
transformation of algae, the method comprising: 

providing at least one alga, immersed in, or adjacent to a 
surface of, a liquid; 

40 providing a mobile system that floats on and is submersible 
in the liquid to a controllable depth, the mobile mecha- 
nism comprising at least first and second environmental 
sensors that sense distinct first and second environmen- 
tal parameter values, pv(t;l) and pv(t;2), respectively, 
45 where at least one of the first and second environmental 
sensors senses at least one of liquid temperature at a 
selected depth within the liquid, temperature of gas adja- 
cent to and above the liquid, liquid pH, liquid salinity, 
liquid turbidity, 0 2 dissolved in the liquid, C0 2 concen- 
50 tration in the liquid, oxidization and reduction potential 
of the liquid, nitrate concentration in the liquid, ammo- 
nium concentration in the liquid, bicarbonate concentra- 
tion in the liquid, phosphate concentration in the liquid, 
light intensity where light enters the liquid, electrical 
55 conductivity of the liquid, and an alga parameter that is 
a measure of growth stage of the alga, and another of the 
at least two environmental sensors senses light intensity 
at a selected depth in the liquid; 
sensing the first environmental parameter value pv(t;l) at 
60 spaced apart first and third sensing times, t=tl and t=t3; 
sensing the second environmental parameter value pv(t;2) 
at a second sensing time, t=t2, that is intermediate 
between the first and third sensing times; and 
transmitting the environmental parameter values, pv(t=tl; 
65 1) and pv(t=t3;l), sensed at the first and third sensing 

times, and the environmental parameter value pv(t=t2; 
2), sensed at the second sensing time, to a data receiver. 
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2. The method of claim 1, further comprising providing a 
locomotion mechanism that is configured to move said 
mobile system along a selected trajectory in said liquid. 

3. The method of claim 1, further comprising providing a 
location sensing mechanism that provides at least first and 
second sets of location observation coordinates of said mobile 
system for at least first and second location observation times, 
respectively, and transmits the at least first and second sets of 
location observation coordinates and the at least first and 
second location observation times to said data receiver. 

4. The method of claim 3, wherein said location sensing 
mechanism provides said at least first and second sets of said 
location observation coordinates of said mobile system by a 
procedure comprising: 

transmitting timed electromagnetic signals from at least 
first, second and third electromagnetic signal sources, 
spaced apart from each other and from said mobile sys- 
tem, at determinable first, second and third location sig- 
nal emission times, respectively, to a location signal 
receiver on said mobile system; and 

estimating at least first and second location coordinates 
(x,y) for a first mobile system location, determined with 
reference to the location signal receiver at a location 
signal receipt time, determined with reference to at least 
one of the first, second and third location signal emission 
times. 

5. The method of claim 3, wherein said location sensing 
mechanism provides at least one of two angular orientation 
coordinates, 0 and <j>, for said mobile system by a procedure 
comprising: 

transmitting timed electromagnetic signals from at least 
first, second and third spaced apart electromagnetic sig- 
nal sources at determinable first, second and third loca- 
tion signal emission times at a first and second, spaced 
apart location signal receivers located on said mobile 
system; 

estimating at least first and second location coordinates 
(x k ,y k ) for each of first and second spaced apart loca- 
tions, corresponding to k==l and k=2, respectively, on 
said mobile system, determined with reference to the 
first location signal receiver and the second location 
signal receiver, respectively, at a location signal receipt 
time, determined with reference to at least one of the 
first, second and third location signal emission times; 
and 

determining at least one of the angular orientation coordi- 
nates, 0 and <|), for said mobile system with reference to 
an angular orientation of a line segment extending 
between first and second locations having the respective 
location coordinates (x^yj and (x 2 ,y 2 ). 

6. The method of claim 3, further comprising: 

providing a locomotion mechanism that is configured to 

move said mobile system along a selected mobile system 
trajectory in said liquid; 

representing the mobile system trajectory as a sequence of 
trajectory location coordinates representing at least first 
and second spaced apart locations, through which the 
mobile system trajectory passes, and corresponding first 
and second location observation times; 

receiving the location observation coordinates at a location 
signal receiver and providing observation location coor- 
dinates for the at least first and second location observa- 
tion times; 

comparing the location observation coordinates with the 
mobile system trajectory location coordinates for the at 
least first and second locations through which the mobile 
system trajectory passes; and 
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when a distance between the location observation coordi- 
nates and the mobile system trajectory location coordi- 
nates, for at least one of the first and second locations 
through which the mobile system trajectory passes, is 

5 greater than a selected threshold distance, issuing a cor- 
rection signal to vary at least one location observation 
coordinate of the mobile system relative to the mobile 
system trajectory. 

7. The method of claim 1, wherein said mobile system is 

1 0 permitted to move in said liquid according to prevailing liquid 

currents in said liquid. 

8 . The method of claim 1 , further comprising: 
providing an elapsed time difference value between said 

15 first sensing time and said second sensing time for said 
first sensed parameter value; 

providing an estimation of a time rate of change of said first 
sensed parameter value; and 

adjusting the elapsed time difference so that this difference 

20 decreases monotonically with increase in the time rate of 

change of said first sensed parameter value. 

9. Hie method of claim 8, further comprising adjusting said 
elapsed time difference so that a change in said first sensed 
parameter value, between at least two consecutive sensing 

25 times, does not exceed a selected fraction f of an average 
value of said first sensed parameter value, where f lies in a 
range 0<f<l . 

10 . The method of claim 9 , further comprising selecting 
said fraction f in a range of 0.01 to 0.20. 

30 1 1 . The method of claim 1 , wherein said first sensed param- 

eter value pv(t;l) and said second sensed parameter value 
pv(t;2) are transmitted with respective first and second 
parameter value reporting frequencies that are different from 
each other. 

35 12. The method of claim 1, wherein said first sensed param- 

eter value pv(t;l) and said second sensed parameter value 
pv(t;2) are transmitted with respective first and second 
parameter value reporting frequencies that are the same. 

13 . The method of claim 1 . further comprising: 

40 receiving at least 2K+1 first sensed parameter values in a 
selected time interval, where K is a selected positive 
integer; 

forming a weighted average WA1 of the at least 2K+1 
received first sensed parameter values; and 

45 replacing each of the at least 2K+1 first sensed parameter 
values by the weighted average AW1. 

14 . The method of claim 1 . further comprising: 
receiving a sequence of at least 2K+1 first sensed param- 
eter values in a selected time interval, where K is a 

50 selected positive integer; and 

selecting a first sensed parameter value that is a median of 
the sequence. 

15 . The method of claim 1 , further comprising: 
receiving at least 2K+1 first sensed parameter values in a 

55 selected time interval, where K is a selected positive 
integer; 

arranging the at least 2K+ 1 first sensed parameter values in 
a sequence that is monotonically increasing or mono- 
tonically decreasing; and 

60 choosing a first sensed parameter value that represents at 
least one of a majority view choice and a plurality view 
choice among the sequence. 

16 . The method of claim 3 , further comprising providing a 
submersion mechanism, attached to said mobile system, that 

65 permits said mobile system to be submerged or raised to a 
selected and continuously variable or discretely variable 
depth. 
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17. The method of claim 16, wherein said mobile system is 
submerged within said liquid by a procedure comprising: 

providing at least first and second bags that are inflatable 
with a selected gas and that are attached at spaced apart 
locations on said mobile system, where the first and 
second bags contain first and second selected amounts, 
respectively, of first and second selected gases; 
when said mobile system is to be raised to said selected 
depth within said liquid, increasing the selected amount 
of the selected gas in at least one of the first and second 
bags; and 

when said mobile system is to be lowered to said selected 
depth within said liquid, decreasing the selected amount 
of the selected gas in at least one of the first and second 
bags. 

18. The method of claim 1, further comprising providing 
said mobile system with an external covering having at least 
one aperture through which an adjacent portion of said liquid 
is exposed to sense at least one of said first and second 
environmental parameters so that at least one of said environ- 
mental parameter values can be measured in said liquid. 

19. The method of claim 18, wherein at least two of said at 
least first and second environmental parameters is sensed 
through said at least one aperture. 

20. The method of claim 18, further comprising providing 
a signal transmitting mechanism, having a signal transmis- 
sion antenna, on an upper portion of said external covering. 

21. The method of claim 18, further comprising providing 
said external covering with a lower portion that is shaped as a 
portion of a prolate spheroid that is partly or wholly sub- 
merged in said liquid. 

22. The method of claim 18, further comprising providing 
said external covering with a lower portion that is shaped as 
an M-hedral surface (M^3) that is partly or wholly sub- 
merged in said liquid. 

23. The method of claim 1, further comprising: 
comparing said sensed first environmental parameter value 

pv(tl;l) for a selected alga, sensed at said first sensing 
time, with a range of acceptable values pvl (ref) of said 
first parameter value; 

when said sensed first parameter value pv(tl;l) is not 
within the acceptable range of values pvl (ref), interpret- 
ing this condition as indicating that growth of the 
selected alga should be discontinued. 

24. The method of claim 1, further comprising: 
submerging said mobile system to a selected positive depth 

in said liquid, corresponding to a selected alga having 
selected location observation coordinates, at spaced 
apart location observation times, t=t4 and t=t5; 
estimating said light intensity in said liquid at the selected 
location observation coordinates at the location obser- 
vation times, t=t4 and t=t5; and 
estimating an amount of growth of the selected alga for at 
least one of the times t4 and t5. 

25. A method of monitoring an environment for growth or 
transformation of algae, the method comprising: 

providing at least one alga immersed in, or adjacent to a 
surface of, a liquid; 

providing N mobile systems (Ng2), numbered n=l , . . . , N, 
where each mobile system floats and/or is submersible 
in the liquid to a controllable depth, and each mobile 
system, number n, comprises at least first and second 
environmental sensors that sense distinct at least first 
and second environmental parameter values, pv(t;l ) and 
pv(t;2), respectively, in common, where at least one of 
the first and second environmental sensors senses at least 
one of liquid temperature, temperature of gas in the 
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liquid, liquid pH, liquid salinity, liquid turbidity, 0 2 
dissolved in the liquid, C0 2 contained in the liquid, 
oxidization and reduction potential of the liquid, nitrate 
concentration in the liquid, ammonium concentration in 
5 the liquid, bicarbonate concentration in the liquid, phos- 
phate concentration in the liquid, light intensity where 
light enters the liquid, electrical conductivity of the liq- 
uid, and an alga parameter that is a measure of growth 
stage of the alga, and another of the at least two envi- 
10 ronmental sensors senses light intensity at a selected 
depth in the liquid; 

sensing the first parameter value pv(t;l;n) at spaced apart 
first and third sensing times, t=tl and t=t3, for each of the 
N mobile systems to provide sensed first parameter val- 
15 ues, pv(tl;l;n) and pv(t3;2;n) for the mobile system 
number n; 

sensing the second parameter value pv(t;2;n) at a second 
sensing time, t=t2, which is intermediate between the 
first and third sensing times, to provide a sensed second 
20 parameter value pv(t2;2;n) for the mobile system num- 
ber n; 

transmitting the values, pv(tl;l;n) and pv(t2;2;n), to a data 
receiver that is spaced apart from at least one of the N 
mobile systems. 

25 26. The method of claim 25, further comprising: 

allowing each of said N mobile systems to move with liquid 
currents in said liquid; and 

providing a location sensing mechanism for each of said N 
mobile systems that senses and provides at least first and 
30 second sets of location observation coordinates (x„(t),y„ 
(t),z„(t), (n=l, . . . , N) for said mobile system number n 
for at least first and second spaced apart location obser- 
vation times, t=tl and t=t2, respectively, and transmits 
the at least first and second sets of location observation 
35 coordinates, and the at least first and second location 
observation times and said observed environmental 
parameter values pv(t=tk;n) (k=l ,2) to a data receiver. 

27. The method of claim 26, further comprising providing, 
at said location data receiver a computer that is programmed: 
40 to receive said at least first and second sets of location 
observation coordinates and said at least first and second 
location observation times, t=tl and t=t2; and 

to estimate, for at least one of said location observation 
times, tl and t2, and for at least one of said environmen- 
45 tal parameter values pvl(t;k;n) (k=l, 2; n=l, . . . , N), a 
continuously variable environmental parameter field of 
values. pv(x,y,z;t;k;field), whose value for said location 
observation coordinates for each of said mobile system 
number n (n=l N), for at least one of said location 
50 observation times, tl and t2, is equal to a corresponding 
environmental parameter value pvl(t;k;n) at the corre- 
sponding location observation coordinates (x„(t),y„(t), 
z„(t)) for at least one of said location observation times, 
tl and t2. 

55 28. The method of claim 27, wherein said computer is 

further programmed: 

to compare said environmental fields of parameter values 
pv(x,y,z;t; field) for t=tl and for t=t2 for at least one 
common set of location coordinates, (x=x c ,y=y £ ,,z=z c ,); 
60 and 

to estimate a rate of change of said environmental param- 
eter value pv(x,y,z;t; field) at the common location coor- 
dinates (x=x c ,y=y c ,,z=z < ,) in a time interval, defined by 

tlStSt2. 

65 29. The method of claim 25, further comprising providing 

each of at least two of said N mobile systems, numbers n=nl 
and n=n2, with a mutual repulsion mechanism so that said 
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systems n=nl and n=n2 repel each other when a distance 
between the at least two systems is less than a selected thresh- 
old distance. 

30. The method of claim 29, further comprising: 
providing, as said mutual repulsion mechanism for said 
mobile system number n=nl, a first magnetic pole 
mechanism having a first selected number Ml of mag- 
netic north pole sources (M1S2), spaced apart on an 
exposed surface of said mobile system number n=nl, 
and having a magnetic south pole source, positioned at a 
location on said mobile system number n=nl that is 
spaced apart from each of the Ml north magnetic pole 
sources; and 

providing, as said mutual repulsion mechanism for said 
mobile system number n=n2, a second magnetic pole 
mechanism having a second selected number M2 of 
magnetic north pole sources (M2 = 2 ), spaced apart on an 
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exposed surface of said mobile system number n=n2, 
and having a magnetic south pole source, positioned at a 
location on said mobile system number n=n2 that is 
spaced apart from each of the M2 north magnetic pole 
5 sources, 

whereby at least one of the Ml magnetic north pole sources 
on the mobile system number n=nl and at least one of 
the M2 magnetic north pole sources on the mobile sys- 
tem number n=n2 will repel each other, when the at least 
to one of the Ml magnetic north pole sources is within said 
threshold positive threshold distance of the at least one 
of the M2 magnetic north pole sources. 

3 1 . The method of claim 24, further comprising using PAM 
fluorometry to estimate said light intensity in said liquid at 
15 said selected location observation coordinates. 



